Doping of graphene adsorbed on the a-SiC>2 surface 
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We have performed an ab initio theoretical investigation of graphene sheet adsorbed on amorphous 
SiC>2 surface (G/a-Si02). We find that graphene adsorbs on the a-SiC>2 surface through van der 
Waals interactions. The inhomogeneous topology of the a-SiCh clean surface promotes a total charge 
density displacement on the adsorbed graphene sheet, giving rise to electron-rich as well as hole-rich 
regions on the graphene. Such anisotropic distribution of the charge density may contribute to the 
reduction of the electronic mobility in G/a-SiC>2 systems. Furthermore, the adsorbed graphene sheet 
exhibits a net total charge density gain. In this case, the graphene sheet becomes n-type doped, 
however, with no formation of chemical bonds at the graphene-SiCh interface. The electronic charge 
transfer from a-SiC>2 to the graphene sheet occurs upon the formation of a partially occupied level 
lying above the Dirac point. We find that such partially occupied level comes from the three-fold 
coordinated oxygen atoms in the a-SiC>2 substrate. 
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Considerable progress has been made addressing the 
electronic and structural properties of graphene, how- 
ever, the development of nanodevices based on graphene 
is still in its infancy. Graphene is a semiconductor ma- 
terial with zero energy bandgap. The valence and the 
conduction bands are connected through a linear energy- 
momentum relation at the K points lying at the corners of 
the hexagonal Brillouin zone (Dirac points) QJ . The very 
high electronic mobility, up to 200,000 cm 2 /(V.s), makes 
graphene a quite suitable material to the development 
of new electronic (nano) devices. Indeed, there are ex- 
perimental realizations of field effect transistors (FETs) 
based on graphene sheets In those systems, SiC>2 

has been used as the gate dielectric material. Lemme 
et al. [2] verified a reduction of the electron mobility 
upon the graphene interaction with the Si02 gate, while 
Romero et al. [f| obtained a n-type doping for graphene 
in contact with a SiC>2 gate. In this case, the authors 
attributed the presence of surface states to the electronic 
charge transfer from SiO<2 to the graphene sheet. In con- 
trast, there are experimental results indicating a p-type 
doping for graphene on SiC>2 gate @, § ■ Very recently n- 
and p-type doping of graphene, lying on SiC>2 substrate, 
has been tuned through a suitable graphene-SiC-2 "inter- 
face engineering" Q, where foreign molecules have been 
placed between the graphene sheet and the SiC>2 surface. 
Indeed, there are experimental and theoretical studies 
addressing "interface engineering" applied to modify the 
electronic and structural properties of carbon nanotubes 
adsorbed on Si02 [UJ as well as on the other semicon- 
ductor surfaces like Si(001) [HI and InAs(lll) Q3. 



The most of the experimental results indicate that 
there are no covalent bonds between the graphene sheet 
and the amorphous Si02 substrate (a-SiO^). It has 
been suggested a physisorption process [l3l - fl5| where van 
der Waals (vdW) interactions play an important role. 
Graphene sheet feels the corrugations of the a-SiC-2 sur- 
face, giving rise to electron-rich as well as hole-rich re- 
gions on graphene, so called "electron- hole puddles" [llij . 
Based upon the aforementioned findings, the understand- 
ing of the electronic and the energetic properties, within 
an atomic scale, of graphene on a-Si02 surface is a im- 
portant issue to the examined by using ab initio total 
energy calculations. 

In this letter we present ab initio calculations of 
graphene adsorbed on the a-SiC>2 surface, G/a-Si02- The 
calculations were performed within the density functional 
theory (DFT) including the long-range vdW interactions. 
Our total energy results reveal that vdW interactions rule 
the graphene adsorption on the a-Si02 surface. The inho- 
mogeneous topology of a-Si02 gives rise to electron-rich 
and hole-rich regions on the adsorbed graphene sheet, in 
accordance with the experimentally observed "electron- 
hole puddles" for graphene lying on a-Si02. However, 
there is a net charge density gain on the graphene sheet. 
The electronic charge transfer from a-Si02 to graphene 
occurs due to the formation of a partially occupied state 
above the Dirac point. We find that such partially oc- 
cupied state is attributed to the formation of three-fold 
coordinated O atoms within a-Si02- 

The amorphous structure was generated through ab 
initio molecular dynamics (MD) simulations based on the 
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FIG. 1: Atomic structure of graphene adsorbed a-Si02 sur- 
face, G/a-Si02 (a) side-view, (b) top-view, details of the equi- 
librium geometry of O t hreo-foid (c). In (a) the bottommost O 
dangling bond was saturated by a H atom and the dashed re- 
gion indicate a periodic unit cell. In (c) the atomic distances 
are in A. Filled (Empty) circles represent the O (Si) atoms, 
(d) Total charge density fluctuation (Ap) within a region of 
±0.5 A with respect to the (equilibrium) vertical position of 
the n-type doped graphene sheet. Within the parentheses we 
present Ap for the neutral G/a-Si02 graphene sheet. Ap is 
in unit of 10 14 e/cm 2 . 



DFT approach as implemented in the VASP code [Pt1 - [19| . 
Initially we obtained the a-Si02 bulk phase. We have 
used a supercell with 96 atoms and we made a suitable 
choice of the lattice constants of the supercell, 12.47 and 
12.96 A along the x and y directions, respectively, in or- 
der to match with the calculated equilibrium lattice con- 
stant of an isolated graphene sheet (Fig.[T]). While the su- 
percell volume was minimized along the z direction, thus, 
avoiding any artificial strain in the a-SiC>2 bulk. We next 
include a vacuum region of —12 A parallel to the z direc- 
tion, within the slab method, and then we repeat the MD 
simulations in order to get an a-SiC>2 clean surface. In 
Ref. [20| we present details on the generation procedure of 
amorphous SiC>2 bulk structure. Once we obtained a well 
described a-SiC>2 surface, we start to investigate the equi- 
librium geometry, energetic stability, and the electronic 
properties of a graphene sheet adsorbed onto a-SiC>2 sur- 
face. Here the calculations were performed by using the 
DFT approach as implemented in the SIESTA code [Il[ . 
The generalized gradient approximation due to Perdew, 
Burke and Ernzerhof (PBE) [22j was used, and the van 
der Waals interaction was described within a semiempir- 
ical approach, following the Grimme formula 23]. The 
electron-ion interactions were calculated by using norm- 



conserving pseudopotentials [24]. All the atomic posi- 
tions of the G/a-Si0 2 system were relaxed by using the 
conjugated gradient scheme, within a force convergence 
criterion of 20 me v/A. 

Figure QJa) presents the structural model for the ph- 
ysisorbed G/a-Si0 2 system, where we have considered 
nine G/a-Si02 configurations on two different a-Si02 
surface structures. At the equilibrium geometry we 
find a graphene - a-Si02 equilibrium vertical distance 
(cfc-SiOs) °f 3.64 (3.30) A for an a-SiC>2 surface cor- 
rugation of 1.7 (1.5) A [11]. Although the relatively 
short distance between the graphene sheet and the top- 
most O atoms, we do not find the formation of C-0 
chemical bonds. The strength of graphene adsorption 
on a-Si02 can be measured by comparing the total en- 
ergies of the isolated systems, viz.: graphene sheet and 
a-Si02 surface, with the total energy of the graphene 
adsorbed system, G/a-Si02- We find that the forma- 
tion of G/a-Si02 is an exothermic process with an ad- 
sorption energy (E ads ) of 6.3±0.4 meV/A 2 Our adsorp- 
tion energy and equilibrium geometry results are in good 
agreement with the experimental estimative of Ishigami 
et ai, E ads = 6 meV/A 2 , and d G -Si0 2 = 4.2 A Q3- Re- 
cent theoretical calculations, within the DFT-LDA ap- 
proach with no inclusion of vdW interations, obtained 
E ads _ 2 me y/A 2 and G?G-Si0 2 = 3.6 A for graphene 
on a-Si02 @. In contrast, other experimental measure- 
ments obtained dc-Si0 2 °f 5.8 A for G/a-Si02, which 
may indicate the presence of foreign elements at the G- 
Si02 interface We next calculate E ads by turning off 
the vdW contribution from our total energy results. In 
this case, the formation of G/a-Si02 becomes no longer 
exothermic, E ads = —0.4 meV/A 2 . Thus, showing that 
graphene sheet adsorbs on a-Si02 is mediated by vdW 
interactions. 

Recent experimental findings indicate a corrugation 
of —5 A on the a-Si02 surface, for a lateral scale of 
-100 A [r|- Meanwhile in Ref. [3 the authors sug- 
gest that the corrugation of the adsorbed graphene sheet 
should correspond to around 60% of the one of the a-Si02 
surface. They measured a graphene corrugation of 1.9 A 
and estimate a deformation energy of 1 me V/A 2 . We ob- 
tained a very small corrugation (up to 0.3 A) and defor- 
mation energy (0.7 meV/A 2 ) for the adsorbed graphene 
sheet, whereas our simulated a-Si02 surface exhibits a 
corrugation of —1.6 A. Here we would be able to describe 
the corrugation of G/a-Si02 in a suitable way by increas- 
ing the size of the surface unit cell, and performing an 
additional total energy minimization of the surface area. 
In this case, we may find higher values of E ads , since the 
graphene-Si02 contact area will increase, however, such 
increase of E ads will be limited by the induced strain 
energy upon the deformation of the graphene sheet. 

It is worth to note that graphene adsorbed on a-Si02 
exhibits a quite different picture when compared with 
graphene on a-quartz. In the latter, due to the presence 
of oxygen dangling bonds, there is the formation of C- 
O chemical bonds (with E ads = 300 meV/A 2 ) on the 
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O-terminated surface [2g|. In addition, it has been sug- 
gested a p-type doping of the adsorbed graphene sheet. 
Whereas in a-SiC>2, which also exhibits an O-rich surface 
[Fig. EI a)], we do not find the formation C-0 chemical 
bonds. In our simulated annealing process all the sur- 
face O atoms become twofold coordinated, suppressing 
the (energetically unfavorable) O dangling bonds. Very 
recent DFT-LDA calculations indicate that the graphene 
interaction with the (0001) surface of a-quartz is weak, 
with no C-0 chemical bonds [271 ]. 

Figure [ljd) depicts the net charge density displace- 
ment (Ap) on the adsorbed graphene sheet. Ap can be 
written as, 

Ap = p[G Si o 2 ] - p[G]. 

Where p[Gsio 2 ] an d p[G] represent the planar averaged 
total charge density of the adsorbed (G/a-SiO^) and iso- 
lated graphene sheets, respectively. The planar average 
was performed by considering a region of 0.5 A above and 
below the graphene sheet. The last term, p[G], was cal- 
culated keeping the same equilibrium geometry as that 
of G/a-SiC-2. We find graphene regions with Ap > 
as well as Ap < 0, electron- and hole-rich regions, re- 
spectively. Thus, Fig. QJd) indicates that there is a total 
charge density fluctuation on the graphene sheet, being 
in accordance with the recent experimental verification 
of "electron-hole puddles" in G/a-Si02 systems flrl l28l]. 
Here we can infer that the inhomogeneous a-Si02 surface 
topology rules the total charge density fluctuation of the 
adsorbed graphene sheet. Such fluctuation on the elec- 
tronic charge density distribution on the graphene sheet 
should play an important role on the experimentally ver- 
ified reduction of the electronic mobility for graphene ly- 
ing on a-Si02 surface 0). 
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FIG. 2: Electronic band structure (a) and the projected den- 
sity of states (b) of G/SiC>2 system. Dark shaded regions in- 
dicate the total density of states, and the light shaded region 
indicate the projected density of states of a-Si02 surface, and 
the solid line indicates the projected density of state of (ad- 
sorbed) graphene sheet. The inset in (a) present the energy 
dispersion of the partially occupied level cf. 

There is a net electronic charge density gain of about 



3 xlO 13 e/cm 2 , Fig.QJd), thus, indicating a n-type dop- 
ing on the adsorbed graphene sheet. Our electronic 
band structure calculation, Fig.EJa), confirms statement 
above. The shape of the Dirac cone has been maintained, 
since the graphene sheet weakly interacts with the a-Si02 
substrate. Due to the electronic charge transfer from a- 
Si02 to the graphene sheet, the Dirac point lies below 
the Fermi energy (Ep). The n-type doping occurs due 
to the presence of a partially occupied state (cf) below 
the conduction band minimum (CBM) pinning the Fermi 
level above the Dirac point, see inset of Fig. [2][a). Fig- 
ure[2{b) presents the density of states (DOS) of G/a-Si02 
(gray region), and separately the contribution from the 
adsorbed graphene sheet (solid line), and the a-Si02 sur- 
face (light gray region). We find, in contrast with previ- 
ous investigations [6J, that cl comes from the formation 
of three-fold coordinated O atoms (Othrcc-foid) in a-Si02- 
Those Othrcc-foid sites may occur when we have a slightly 
hight local concentration of Si atoms or nearby oxygen 
vacancies in a-Si02 [29-32]. The Si-O t hrcc-foid bond 
lengths (1.82-1.86 A) are stretched in comparison with 
the ones of the two fold coordinated O atoms, ^f .66 A, 
Fig. []Tc) . Those larger bond lengths reduce the binding 
energy of the electronic states along Si-Othree-foid, giv- 
ing rise to the partially occupied state, cl, just below the 
CBM and pinning the Fermi level. We find that the elec- 
tronic states of cf are localized on the Si atoms nearby 
Othree-foid- In order to verify the adequacy of the pro- 
posed origin of cl, we calculate the electronic structure 
of G/a-Si02 by suppressing the presence of Othrcc-foid, *■ 
e. we have only two fold coordinated O atoms within the 
a-Si02 slab. In this case, we do not find the formation 
of cf , and the Dirac point lies on the Fermi level. That 
is, the adsorbed graphene sheet becomes neutral. Those 
results support the proposed role of Othrcc-foid on the 
non-covalent n-type doping of graphene on a-Si02 sub- 
strate. It is worth to note that, for the neutral graphene 
sheet, we find the same picture for the electron- and hole- 
rich regions, as depicted in Fig.[ljd), however with differ- 
ent values of Ap. The total charge density fluctuation on 
the neutral graphene is about 25% smaller in comparison 
with the (n-type) doped ones. 

Finally, we have considered the presence of oxygen va- 
cancies (Ov) in a-Si02- We have examined three plausi- 
ble Ov configurations, on the topmost surface site, and 
on the subsurface sites, indicated as V1-V3 in Fig. [TJa), 
respectively. We verify that those vacancies do not give 
rise to a donor level on the a-Si02 surface. For VI and V2 
we find the formation of Si-dimers on the surface (Si-Si 
bond length of 2.24 A) inducing the formation of occu- 
pied (empty) electronic states near the valence (conduc- 
tion) band maximum (minimum) . Whereas V3 gives rise 
to a deep occupied state within the a-Si02 bandgap. 

In conclusion, we find that vdW interactions rules the 
graphene adsorption on the a-Si02 surface. Through a 
charge density map on the adsorbed graphene sheet we 
verify that the inhomogeneous a-Si02 surface topology 
promotes an electronic charge density displacement on 
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the graphene, giving rise to the experimentally verified 
"electron- hole puddles" 0, [H[ . Such electron- rich and 
hole-rich regions should contribute to the experimentally 
verified reduction on the electronic mobility of graphene 
sheets lying on a-SiC>2 surfaces. We find a net charge den- 
sity gain on the adsorbed graphene sheet, thus, character- 
izing a n-type non-covalent doping on the graphene. The 
electronic charge transfer from a-SiC>2 to the graphene 
sheet comes from the formation of a partially occupied 
state above the Dirac point, which is attributed to the 
formation of three-fold coordinated oxygen atoms in a- 
SiC>2. Those three-fold coordinated oxygen atoms may 



occurs upon slightly hight local concentration of Si atoms 
or nearby oxygen vacancies. 



Acknowledgments 

The authors thanks Dr. Renato Pontes for fruitful dis- 
cussions, the financial support from the Brazilian agen- 
cies CNPq, CAPES, FAPEMIG, FAPESP, CNPq/INCT 
and the computational support from CENAPAD/SP. 



[1] A. C. Neto, F. Guinea, A. H. C. Neto, N. M. R. Peres, 

K. S. Novoselov, K. Novoselov, A. K. Geim, and A. Geim, 

Rev. Mo.d Phys. 81, 109 (2009). 
[2] M. C. Lemme, T. J. Echtermeyer, M. Bauss, and H. Kurz, 

IEEE Elect. Dev. Lett. 28, 282 (2007). 
[3] X. Li, X. Wang, L. Zhang, S. Lee, and H. Dai, Science 

319, 1229 (2008). 
[4] X. Wang, Y. Ouyang, X. Li, H. Wang, J. Guo, and 

H. Dai, Phys. Rev. Lett. 100, 206803 (2008). 
[5] J. B. Oostinga, H. B. Heersche, X. Liu, A. F. Morpurgo, 

and L. M. K. Vandersypen, Nature Materials 7, 151 

(2008). 

[6] H. Romero, N. Shen, P. Joshi, H. R. Gutierrez, S. A. 
Tadigadapa, J. O. Sofo, and P. C. Eklund, ACS Nano 2, 
2037 (2008). 

[7] S. Berciaud, S. Ryu, L. E. Brus, and T. F. Heiz, Nano 

Lett. 9, 346 (2009). 
[8] S. M. Song and B. J. Cho, Nanotechnol. 21, 335706 

(2010). 

[9] R. Wang, S. Wang, D. Zhang, Z. Li, Z. Fang, and X. Qiu, 
ACS Nano 5, 408 (2011). 
[10] J. S. Soares, A. P. M. Barboza, P. T. Araujo, N. M. B. 
Neto, D. Nakabayashi, N. Shadmi, T. S. Yarden, A. Is- 
mach, N. Geblinger, E. Joselevich, et al., Nano Lett. 10, 
5043 (2010). 

[11] W. Orellana, R. H. Miwa, and A. Fazzio, Phys. Rev. Lett. 

91, 166802 (2003). 
[12] Y. H. Kim, M. J. Heben, and S. B. Zhang, Phys. Rev. 

Lett. 92, 176192 (2004). 
[13] E. Stolyarova, K. T. Rim, S. Ryu, J. Maultzsch, P. Kim, 

L. E. Brus, T. F. Heinz, M. S. Hybertsen, and G. W. 

Flynn, PNAS 104, 9209 (2007). 
[14] M. Ishigami, J. H. Chen, W. G. Cullen, M. S. Fuhrer, 

and E. D. Williams, Nano Lett. 7, 1643 (2007). 
[15] A. Sinitskii, D. V.Kosynkin, A. Dimiev, and T. J. M, 

ACS Nano 4, 3095 (2010). 
[16] J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. 



Smet, and K. Von Klitzing, Nat. Phys. 4, 144 (2008). 
[17] G. G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 
[18] G. G. Kresse and J. Hafner, J. Phys. 6, 8245 (1994). 
[19] G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 

(1996). 

[20] W. L. Scopel, A. J. R. da Silva, and A. Fazzio, Phys. 
Rev. B 77, 172101 (2008). 

[21] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon, and D. Sanchez-Portal, J. Phys.: Con- 
dens. Matter 14, 2745 (2002). 

[22] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[23] S. Grimme, J. Comput. Chem. 27, 1787 (2006). 

[24] N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

[25] The a-Si02 surface corrugations was inferred by the cal- 
culation of the variance of vertical positions of the top- 
most O atoms, while for the adsorbed graphene sheet, 
we considered the variance of vertical positions of the 
graphene C atoms. 

[26] Y. J. Kang, J. Kang, and K. J. Chang, Phys. Rev. B 78, 
115404 (2008). 

[27] T. C. Nguyen, M. Otani, and S. Okada, Phys. Rev. Lett. 

106, 106801 (2011). 
[28] C. Stampfer, F. Molitor, D. Graf, K. Ensslin, A. Junger, 

C. Hierold, and L. Wirtz, Appl. Phys. Lett 91, 241907 

(2007). 

[29] Z. Y. Lu, C. J. Nicklaw, D. M. Fleetwood, R. D. 
Schrimpf, and S. T. Pantelides, Phys. Rev. Lett. 89, 
285505 (2002). 

[30] V. O. Sokolov and V. B. Sulimov, J. NosnCryst. Sol. 217, 
167 (1997). 

[31] B. P. Feuston and S. H. Garofalini, J. Chem. Phys. 91, 
564 (1989). 

[32] K. Winer and L. Ley, Phys. Rev. B 37, 8363 (1988). 



